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ABSTRACT 

The central densities of dark matter (DM) halos are much lower than predicted in cold DM models 
of structure formation. Confirmation that they have cores with a finite central density would allow us 
to rule out many popular types of coUisionless particle as candidates for DM. Any model that leads to 
cusped halos (such as cold DM) is already facing serious difficulties on small scales and hot DM models 
have been excluded. Here I show that fermionic warm DM is inconsistent with the wide range of phase 
space densities in the DM halos of well-observed nearby galaxies. 

Subject headings: galaxies: kinematics and dynamics - galaxies: halos - galaxies: dark matter - 
galaxies: formation 



1. INTRODUCTION 

The well-known mass discrepancies in galaxies (Bosma 
1978; Broeils 1992; Verheijen 1997) and in clusters of 
galaxies (Zwicky 1937; Carlberg, Yee & Elhngson 1997; 
Tyson, Kochanski & Dell'Antonio 1998) are usually taken 
to imply the existence of a large fraction of invisible, or 
"dark," matter (DM) in the universe. A popular candidate 
is Cold Dark Matter, for which the DM particles, what- 
ever they are, are almost at rest with respect to the Hubble 
flow in the early universe. CDM is often imagined to be a 
heavy, non-baryonic relic particle from the early universe 
which has essentially only gravitational interactions with 
itself and with normal, or baryonic, matter. 

The CDM model has been studied intensively for 
twenty years and now has many well-worked out predic- 
tions for the formation of structure in the universe {e.g. 
Bertschinger 1998). The broad-brush impression one now 
has is that the currently favored ACDM model boasts 
a considerable degree of success in predicting large-scale 
structure {e.g. Pearce et al. 1999; Bahcall et al. 1999). But 
it has become apparent in recent years that the predictions 
of almost any flavor of CDM are seriously at variance with 
the observed properties of galaxies because the central den- 
sities of collapsed objects and fragments are predicted to 
be too high (§|). 

A number of authors have therefore begun to explore 
variations of the CDM model; most favor a modification 
to the properties of the DM particles rather than the al- 
ternative - a change to the law of gravity. The simplest 
is the warm DM matter model {e.g. Colombi, Dodelson 
& Widrow 1996; Sommer-Larsen & Dolgov 1999; Hogan 
1999) in which streaming of the DM particles in the early 
universe suppresses small-scale power in the fiuctuation 
spectrum. In addition, WDM particles in halos of greater 
volume density must have larger velocity spreads, because 
of Liouville's theorem, thereby precluding strong density 
gradients. 

Simulators of the WDM model suppress small-scale 
power in the fluctuation spectrum but generally ignore 
the initial flnite velocity spread, which is difficult to in- 



clude without wrecking the quiet start. Since the DM in 
these simulations still has infinite phase space density, the 
resulting halo profiles have density cusps resembling those 
which form in CDM {e.g. Moore et al. 1999; Colfn et al. 
2000). 

2. PHASE SPACE DENSITY CONSTRAINT 
2.1. Prediction 

If the DM particle decoupled from thermal equilibrium 
at some early epoch, then its phase space density today 
can be predicted. Fermionic DM would have a flnite ini- 
tial maximum / = /max, which leads the well-known con- 
straint on neutrino masses (Tremaine & Gunn 1979). 

Hogan (1999) suggests that halos today may reflect, /, 
the average over the initial velocity distribution. But /max 
is unbounded for bosons, making a single value of / rather 
poorly deflned (Madsen 2000). Thus phase space density 
constraints are much weaker for bosonic WDM; careful 
simulations are needed to predict the structure of halos in 
this model, but the small fraction of high / material makes 
mild cusps seem likely. 

The flne-grained phase space density, /fine, of a truly 
coUisionless fluid is strictly conserved, but the maximum 
coarse-grain density, /coarse, can decrease during violent 
evolution, such as a collapse or merger. Simulations, which 
by their nature can track only /coarse, have found, however, 
that the maximum /coarse barely changes in even quite 
violent collapses (May & van Albada 1984) or mergers 
(Farouki, Shapiro & Duncan 1983; Barnes 1992). These 
tests were with systems that have flnite /max at the outset, 
whereas /coarse can, and generally does, decrease during 
the collapse or merger of systems with unbounded /max- 
Hernquist, Spergel & Heyl (1993) merged systems with 
cusped density proflles having unbounded /max; they de- 
flned an "average" phase space density over the inner half 
of the mass distribution fi/2, which decreased during ma- 
jor mergers. Dalcanton & Hogan (2000) construct a heuris- 
tic argument that /coarse should decrease during mergers, 
but like Hernquist et al., their argument applies to the av- 
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erage value over the inner part of the halo, which probably 
does decrease, rather than the specific central value which, 
if originally finite, remains almost unchanged. 

Processes such as baryonic cooling and smooth infall, or 
even feedback from star formation, do not change / for the 
DM. Scattering of halo particles by dense clumps of bary- 
onic matter could change / locally, but simulations with 
a baryonic component {e.g. Navarro & Steinmetz 2000) 
do not appear to lead to significantly different halo pro- 
files from those in which baryonic processes are omitted. 
Dynamical friction in barred galaxies also hardly affects 
the structure of the inner halo, even when it is intolerably 
fierce (Debattista & Sellwood 2000). 

Thus /max for halos today should all cluster around a 
common value if the DM particle is a coUisionless fermion 
having some fixed velocity spread in the early universe. 

2.2. Phase Space Densities of DM Halos 

Fermionic WDM halos should today approximate cored 
isothermal spheres having finite central densities. The 1- 
D velocity dispersion in an isotropic, isothermal sphere 
is simply a = 2~^/^14at, where 14iat is the flat circular 
speed from the halo (Binney & Tremaine 1987 §4. 4b). A 
finite central density, po, yields a rotation curve with a 
central slope dV/dR = (47rG/9o/3)^^^. For the popular 
halo fitting function p = pq/{1 + r'^/r^), the asymptotic 
velocity, core radius and central density are related as 
Vflat — (47rG'po)^^^^o- Other functional forms {e.g. Evans 
1993), departures from sphericity, or from velocity isotropy 
merely introduce correction factors of order unity. Thus 
measurements of Vflat and of either po or rg for the halo 
allow us to estimate /max ^ Po/o'^- 

Rotation curves of galaxies include contributions from 
components other than DM, of course. HI data that ex- 
tend well outside the optical part of the galaxy provide 
reasonably firm values for Vflat- Rotation curves that are 
well-resolved in the inner parts provide an upper limit to 
po, which may be close to the actual value if the baryonic 
contribution is small, as in low surface brightness systems 
(LSBs). In many cases, however, either Vflat or po for 
the DM halo depends on the decomposition of the rota- 
tion curve into the separate contributions from luminous 
and dark matter, which is generally controversial. I as- 
sume maximum disk models, and discuss below how this 
assumption affects the values. 

Figure |l| shows estimates of /max ~ Po/c'^ in a number 
of well-observed galaxies. The data and their sources are 
summarized in Table |l] with the exception of the points 
for the Draco and Ursa Minor dwarf spheroidal galaxies. 
Broeils (1992) gives large formal values of Vflat ^ V^ax for 
three galaxies (noted in Table |l|) ; in these cases, I conser- 
vatively determine cr from VJnax instead of Vflat ■ 

The pluses in Figure ^ are for the Draco and Ursa Mi- 
nor dwarf spheroidal galaxies, which both have a stellar 
velocity dispersion 10 km s~^ (Armandroff, Olszewski 
& Pryor 1995). The uppermost point assumes a simple 
mass-foUows-light King model (Binney & Tremaine 1987, 
§4.4) with an estimated King radius for Draco of 150 pc 
(Pryor & Kormendy 1990) that yields a central density of 
0.7 M0 pc~'^ (Pryor, private communication). The DM 
halo could have a larger rg and a than the stars, however, 
implying a lower /max- We place an extreme lower bound 
on /max, which is essentially the same for both Draco and 



UMin., as follows: We adopt the lower bound po = 0.2 
Mq pc""^ (Olszewski 1998), and the argument (Gerhard 
& Spergel 1992) that their masses must be < 10^° Mq. 
Treating this gigantic, low-density halo as a Wq — 1 King 
model, we obtain ctdm ^ 130 km s~^. Slightly more ex- 
treme models could be imagined, but M/Ly ^ 30 000 for 
this model already! A more resonable lower bound might 
be to adopt M/Ly ~ 500, on the high side for galaxy clus- 
ters {e.g. Carlberg et al. 1997; Tyson et al. 1998), which 
would require ctdm ^ 32 km s^^. The values of /max for 
this, and the above more extreme, model are shown by the 
lower pluses in Figure |l|. 

While maximum disk models are not universally ac- 
cepted, most of the points in Figure |l| would not move 
much if this assumption were dropped. Three of the LSBs 
from Swaters, Madore & Trewhella (2000) are plotted 
(squares) for both their "maximum disk" and "no disk" 
fits for which /max differs by at most factor 100. The DM 
halo for a "no disk" fit to NGC 3198, the proto-typical 
galaxy for DM studies, has the same V^flat but a core some 
5-6 times smaller (van Albada et al. 1985), increasing /max 
by the square of this factor. While these are substantial 
uncertainties, they are small in relation to the total range. 

2.3. Implication 

There are three possible conclusions from this Figure. 
First, the range of /max could be mostly due to errors. 
Clearly some invidual points could be in error by a cou- 
ple of orders of magnitude, but others are more certain. 
The spread is admittedly greatly increased by the highest 
point for the dwarf spheroidals, which is not well deter- 
mined because we observe only the stars at the centers of 
the halos. This point can be brought into the main cluster 
in the Figure, but only by adopting a truly extreme model. 
Similarly, the lower points could be moved up by adopting 
sub-maximal disk models. If one wishes to argue that the 
entire spread is due to measurement errors, and that the 
DM halos manifest a characteristic value for /max, that 
value corresponds to a thermal relic fermion with a mass 
of <, 100 eV, i.e. hot DM (HDM), which has been rejected 
previously {e.g. Gerhard & Spergel 1992; Cen & Ostriker 
1992). 

Dalcanton & Hogan (2000), who present a similar dia- 
gram, argue that the measurements are not of /max but 
of some average value over the inner halo. If the den- 
sity profiles of all halos rise to a finite central value with 
a similar functional form, however, their quantity would 
also have a characteristic value in WDM. To be consis- 
tent with WDM, they require the halo density to continue 
to rise steeply inside tq, with high / material at the very 
center contributing little to the average over the volume 
inside rg. Their interpretation therefore therefore applies 
to a DM particle that has initially unbounded /max- 

If the values in Figure 1 are indeed measurements of 
/max, most of the seven orders of magnitude spread is enor- 
mously larger than could be allowed if the WDM particle 
were a coUisionless fermion. 

3. HALOS WITH CUSPS 

Phase space density constraints cannot be applied to 
CDM, or strictly to bosonic WDM despite its velocity dis- 
persion, since /max is unbounded. It is the absence of an 
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upper bound to the phase space density in CDM which 
gives rise to density cusps in the collapsed halos (Moore 
et al. 1998; Klypin et al. 2000). (This is also true of the 
cold component of the less-appealing mixed DM models 
[e.g. Kofman et al. 1996].) An unbounded /max is the root 
cause of the difficulties now besetting CDM: 

(1) The "concentration index" of the halo, which has 
a range of values (Bullock et al. 1999), implies a high 
central DM density, that should increase still further as 
the baryons cool and settle to the center. Low luminosity 
galaxies and LSBs are believed to have the largest frac- 
tions of DM, and therefore halos for which compression by 
baryonic infall is least important. Yet the rotation curves 
of these galaxies (Cote et al. 1997; Swaters et al. 2000) 
rise gently, suggesting a low density core, irrespective of 
the M/L ascribed to the baryonic component. The halo 
profiles in bright HSB galaxies are more model-dependent, 
since the rotation curves in these galaxies generally do rise 
quickly (Rubin, Kenney & Young 1997; Sofue et al. 1999). 
It is now clear, however, that the inner part of the rotation 
curves of these galaxies is dominated by the luminous disk 
and bulge; even in these cases, the DM halo has a large, 
low-density core (Debattista & Sellwood 1998) or very low 
concentration index (Weiner et al. 2000). Thus, there is 
little support for high central densities of DM within the 
halos of any type of galaxy. 

(2) The merging hierarchy causes the cooled baryonic 
fraction to lose angular momentum to the halo, making 
disks that are too small (Navarro & White 1994; Navarro & 
Steinmetz 1997). The predicted angular momentum of the 
disk is an order of magnitude less than that observed. The 
problem is only partially ameliorated (MacLow & Ferrara 
1999; Navarro & Steinmetz 2000) if some process (usu- 
ally described as "feedback from star formation" ) prevents 
most of the gas from cooling until after the galaxy is as- 
sembled. 

(3) Navarro & Steinmetz (2000) describe their failure 
to predict the zero-point of the TuUy-Fisher relation as 
a "fatal problem for the ACDM paradigm." They show 
that no matter what M/L is assumed for the disk, the 
predicted circular speed at a given luminosity is too high 
because the halo density is too high. The TuUy-Fisher pre- 
diction may be even worse, since CDM predicts L 
(Dalcanton, Spergel & Summers 1997; Mo, Mao & White 
1998), whereas Verheijen (1997) stresses that when V is 
interpreted as the circular velocity of the flat part of the 
rotation curve, the true relation is closer to L oc y'. 

(4) The results from high-resolution A''-body simulations 
(Klypin et al. 1999; Moore et al. 1999) have revealed large 
numbers of sub-clumps within large DM halos, many more 
than are observed as satellite galaxies. Whether these frag- 



ments threaten the survival of thin disks in the host galaxy, 
remains to be seen. (WDM models are largely motivated 
to avoid this problem by suppressing small-scale power in 
the initial perturbation spectrum.) 

These difficulties may not yet be fatal to CDM, since a 
better understanding of baryonic processes in galaxy for- 
mation could conceivably alter the predictions. It is un- 
clear what process could loosen the high density clumps 
of DM, however; e.g. Debattista & Sellwood (2000) show 
that an intolerable degree of dynamical friction has a very 
mild effect on the halo density profile. 

The first of these problems, which is shared by all DM 
models that predict cusped halos, is key. The highest qual- 
ity data {e.g. Blais-Ouellette et al. 1999) do appear to show 
that the halo density profile rises gently to a finite central 
value. This result is controversial, however, since others 
{e.g. van den Bosch et al. 1999) argue that the observed ro- 
tation curves could be consistent with mild density cusps 
in the halos. More high-quality data should eventually 
settle the issue. 

4. CONCLUSIONS 

Fermionic WDM halos should have well-defined cores 
with a characteristic /max- This prediction is inconsistent 
with the range of values shown in Figure 1 which rules out 
this particular DM candidate. Halos having finite central 
densities are also inconsistent with the density cusps pre- 
dicted by CDM, and possibly also bosonic WDM. Other 
work has excluded HDM. Thus the existence of finite cen- 
tral density cores in halos would rule out most forms of 
simple, collisionless relic particle as DM candidates. It is 
therefore of great importance to obtain tighter observa- 
tional constraints on the inner density profiles of halos. 

It is suddenly popular to hypothesize extra properties 
for the CDM particle in order to soften the collapsed ha- 
los (Spergel & Steinhardt 1999; Peebles 2000; Goodman 
2000; Hu et al. 2000; Kaplinghat et al. 2000; etc.). The 
simple CDM model would become much less attractive if 
an extra ad hoc property were needed to rescue it. 

The nature of DM is increasingly constrained by the 
observed properties of galaxy halos. In particular, any 
successful theory of galaxy formation will need to account 
for the wide range of phase space densities for DM shown 
in Figure 1. 
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NSF grant AST 96/17088 and NASA LTSA grant NAG 
5-6037. 
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Galaxies Plotted in Figure |i| 
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^Fitted core radius in kpc 

^Fitted asymptotic velocity in km 

"Fitted central DM density in IO^'^Mq pc"^ 

'^Observed maximum circular velocity in km s 

"Adopted Viiat reduced to observed Vmax 
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Fig. 1. — Estimates of /max {~ po/a^) for DM in a number of galaxies plotted against the observed Vmax- The pluses and squares are 
described in the text, the circles are from Broeils (1992) and Verheijen (1997), the cross is for NGC 4123, and the triangle is for NGC 5585. 
Table hi gives the raw data and its source for each galaxy. The units of /max are Mq pc~^ / (km s~^)~'^ and km s~^ for Vmax; the right-hand 
axis snows values of the rest mass of an equivalent thermal relic fermion. 



